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I. INTRODUCTION 



The famous Henon-Heiles potential (for details, see Section 2) was introduced to model the motion of a star within 
a galaxy. But it has much larger connotations. 

Mathematically speaking, the model was constructed by adding two terms of third degree to the potential of a planar 
oscillator. Such a mathematical model also appears by expanding the potential corresponding to an integrable system 
(resulting via some canonical transformations applied to the motion of three particles on a circle under exponentially 
decreasing forces) to third degree terms (Boccaletti and Pucacco 1996 lj]; Anisiu and Pal 1999 Q). 

This problem originates in a question largely discussed in the sixties and approached much earlier: the existence 
of a third isolating integral of the motion (besides the integrals of energy and angular momentum). For terminology, 
significance, and endeavours , see (Henon and Heiles 1964 [3| and the references therein). 

Starting with 1957, G. Contopoulos searched for and found some cases of two potentials when such an integral ( -Z3, 
in the classical notation) exists. Quoting him [Hill' numerical results provide ample evidence of the existence of I3 
in quite general potential fields. As concrete exemplifications from astronomy, but not only: slightly elliptical stellar 
clusters or galaxies, Schmidt's model of the Galaxy, noncentral gravitational field of the Earth; also, systems where 
the unperturbed frequency ratio is rational, or potential fields deprived of a symmetry plane. 

Besides these situations, there exist problems of statistical mechanics @, HQ) celestial mechanics @ and quantum 
mechanics [Io| that join a Henon-Heiles-type model. 

In 1964 M. Henon and C. Heiles Q simplified the initial problem and approached it numerically via Poincare 
sections. Their paper and their numerical experiments had a great echo, benefitting today of more than 600 citations 
in the mathematical, physical and astronomical literature. The model bears their names, and the generalizations are 
called Henon-Heiles-type models. 

As a matter of fact, in many Henon-Heiles-type models (and earlier models) the integral of angular momentum 
does not hold. Such models were also tackled via the tools of the theory of dynamical systems. Only one example at 
hand: the study of the collision and escape dynamics in the associated two-body problem [lj, [l^] • 

Coming back to the Henon-Heiles-type model detailed in Section [H] in 1982, T. Bountis et al. [l3| used Painleve 
analysis to identify all integrable cases. They identified three such situations, ruled by the values of the free parameters 
of the model; out of them case (ii) is the most general (see Section Hlf . 

The present paper intends to point out exact solutions for the class of Henon- Heiles-type systems described by the 
case (ii) above. Section HI1 establishes the basic equations for the problem to be investigated. We start from cylindrical 
coordinates (in the motion reduced to a meridian plane), but then we work with a generalized Henon-Heiles-type 
system in configuration-momentum coordinates. In Section Hill we derive the Lax representation for the corresponding 
system. 

Section IIVI provides exact quasi-periodic solutions for the system under consideration. These orbits are expressed 
in terms of Kleinian hyperelliptic functions. 

In Section [V] we search for elliptic periodic solutions of our system. To this end in view, we resort to the method 
proposed in 1J, EE EH] , which allows the construction of periodic solutions in a straightforward way, applying the 
spectral theory for the Schrodinger equation to elliptic potentials. In this way we get the corresponding solution of 
the system under consideration. 

Section IVT1 tackles a generalized Henon-Heiles-type system with (n + 1) degrees of freedom. We apply the methods 
presented in the previous sections to such a system, and point out the exact solutions in this case. 

The final Section IVlII summarizes the main results obtained in the paper and formulates some conclusions. 



II. BASIC EQUATIONS 

We consider here a model for the dynamics of a point mass (a star) within the gravitational field ruled by the 
potential V g of an axially symmetric galaxy. Let the mass of the star be to, and let the cylindrical coordinates we 
shall use be (r,cj>,z). The axis Oz is the axis of symmetry, z is the distance of the star from the reference plane, 
r := yj x 2 + y 2 is the distance between the star and the axis Oz and <ft '■— arctan^ is the polar angle. Here the 
Oxy-pl&ne is a reference plane (in rectangular coordinates) that defines the coordinate r and the angle <f>. 

The Hamiltonian of the model can be written as 

H{jp r ,p z , P4> ,r,z,<j)) = ^- (p 2 r +p 2 z ) + ^-pl + V g {r,z), 

where p r = mr and p z = mz are the linear momenta in the r and z directions respectively, p^ = mrcj) is the angular 
momentum around the symmetry axis, whereas ' = is the derivative with respect to the time t. 
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We have two integrals of motion, the total energy, E, and the angular momentum, I, respectively: 

E = V g (r,z) + f (r 2 + r 2 <j> 2 + i 2 ), (1) 

I = mr(j). (2) 
With the help of the second integral I, we reduce the dynamics of the star to the meridian plane (r, z): 

.. _ dV{r,z) .. _ dV(r,z) 



dr ' dz 

I 2 

2 12 ■ 1 



(3) 



V(r,z) := V g (r,z) + fr 2 ^ = V g (r, z) + 



2mr 2 



This problem originates in a question discussed in the sixties: the existence of a third isolating integral of the motion 
see the introductory section. Starting with 1957, G. Contopoulos [H found some cases of potentials V(r, z) when 
such integral I3 exists. 

In 1964 M. Henon and C. Heiles [|| simplified the problem, canceling all terms of order > 4 in the potential V(r, z). 
The simplest such potential can be read as 

V(r, z) = Cr 3 + irz 2 + ^ (Ar 2 + Bz 2 ) , 

where A, B and C are free parameters. 

Using Painleve analysis, in 1982 T.Bountis, H. Segur and F.Vivaldi [HI extracted all integrable cases of ((3]), see 
also [H H El: 

(i) A = B, C = I (4) 

(ii) C = 1, A and B are arbitrary (5) 

(hi) 16A = B, C = y . (6) 

In the case (i), the equations of motion decoupled in {z + r), (z — r) coordinates and the general solution can be 
expressed via elliptic functions, see for example (20| and references therein. The question about the effective exact 
solution in the case (iii) is still open [2(| • For the case (ii) effective solutions are obtained in [U [22|, HH, [24[ ■ In this 
paper we shall focus our attention to the Henon-Heiles-type systems related to second case (ii). 

Next we will use 01 and qi instead of r and z. We consider a generalized Henon-Heiles-type system with two degrees 
of freedom [H, [H, [H [H, H| : 

qi + 3qj + \^q\ + a q x - ^ = 0, (7) 
<72 + qi Qi - + ~T92 = 0. (8) 



4qi 4 



Its Hamiltonian is 



where q\, q2,Pi,P2 are the canonical coordinates and momenta and ao, ai, 04 are free constant parameters. Moreover 
Hq is related to the Hamiltonian 



Hh = \ (pi + Pi) + ql + \qiql + ^ + ~ + Sg 2 2 ) , (10) 



through the map 



A 

ft = <?i + - 2-B, 92 = 92, 

a = -2,4+12B, ai = -A 2 + 16AB - 48B 2 . (11) 

The function Hh is the Hamiltonian of a classical integrable Hcnon-Hciles system with the additional term /Sq 2 . 
The corresponding equations are equivalent to the ordinary differential equation for travelling wave solutions of the 
fifth-order flow in the Korteweg - de Vries (KdV) hierarchy (2f| . 
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III. LAX REPRESENTATION 



Next we will derive (2 x 2) matrix Lax representation for the generalized Henon-Heiles system ([9]). The Lax 
representation has the form [30, HH 

L = [M(t, A), L(t, A)], L = ( £ _^ V M = ( ° (12) 



where C7, W, Q are [3J, |3J, [3J see also [35j]: 



tf(t, A) - F(t, A) = A 2 + -^A - 

2 lb 

V(*,A) = -^-F(t, A) = -7P1A+ -U2P2, 
2 4 lb 

W(f, A) = -If + QF = A 3 - + ia )A 2 

+ ( 4 91 + lb 92 - lb" 1 + 8 ao9l)A + lT 2 + 64^|' 



Q(t, A) = A - qi - -a . 
The corresponding algebraic curve is det(L(t, A) — |J) = 



^ 2 = 4A 5 - a A 4 - i ai A 3 + ii/oA 2 + ~#iA - -J-a 4 . (13) 
4 2 8 25b 



It is easy to derive a second integral H\: 



1 2 2 1 4 _ Oo 2 1 fl l 2 

l 4g 2 4 



u 2 ^ 2„2 „4 , 11 >"U u , «u 2 

■ffl = P 2 9l -P1P2Q2 - T;1lQ2 ~ ol2 + fl 4— ~ —91 92 + ~ <72- 



IV. EXACT QUASI-PERIODIC SOLUTIONS 

In this section we g ive the trajectories of the system under consideration in terms of Kleinian hyperelliptic functions 
(see, e.g., [3bl l37l I38LI391] ) . being associated with the real algebraic curve of genus two (|13|). which can be also written 
in the form 



v 2 



= 4[](A-A 4 )=4A 5 + ^a fc A fe , (14) 

2 = fc = 

where \ ^ Ai are branching points and 

a 4 = -a , a 3 = -~<2i, a 2 = ~H , 

(*i = ~H 1 , a = -— a 4 . (15) 
8 25b 

At all real branching points the closed intervals [A2j_i, \2i],i = 1, 2 will be referred further as lacunae pfl.l4l|. Let 
us equip the curve with a homology basis (01, a 2 ; bi, b 2 ) € Hi(K,Z) and fix the basis in the space of holomorphic 
differentials. 

The exact integration of the system {?} and JSJ) reduces to the solution of the Jacobi inversion problem in the 
following form 

X 2 -p 22 {u)X-p 12 {u) = 0, (lb) 
that is, the pair /12) is the pair of roots of p^|) . So we have 

P22{u) = fll +/i2, Pl2(u) = -/J,lfl2- (17) 
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Let us introduce finally the Baker- Akhiezer function, which, in the framework of the formalism developed, is 
expressible in terms of the Kleinian a- function as follows [38| : 

*(A,tt)= ^expiy dr T u , (18) 

where A is arbitrary and u is the Abel image of arbitrary point (y\ 1 fix) x (1/2, /Lt 2 ) € K X K. It is straightforward to 
show by the direct calculation, based on the relations for three and four- index p-functions [38|], that ^ (A, u) satisfies 
the Schrodinger equation 

- 2p 22 (u)) *(A, u) = (x + \a^j *(A, u) (19) 
for all {v,n). The solutions of 0, flSJ have the following form in terms of Kleinian functions p 22 (w), pii{u) [32l l33l. 

MM 

qi = -2p 22 (u), q\ = 16pi 2 (u). 



V. ELLIPTIC PERIODIC SOLUTIONS 



In this section we follow the method proposed in [14j, [15|, lla, 143J , which allows us to construct periodic solutions 
of 0, (H|) in a st raig htforward way based on the application of spectral theory for the Schrodinger equation with 
elliptic potentials (Tj, EE EE 0, H| S3, H EE HE M, H H Hi Hi H . We start with the equation QTJJ) for Baker 
function \I/(A, it). We assume, without loss of generality, that the associated curve has the property 04 = 0. To make 
this assumption applicable to the initial curve of the system (0, (J8| being derived from the Lax representation, we 
undertake the shift of the spectral parameter 

A — ► A - ao/20, (20) 

Consider genus 2 Lame potential u = 6p(t + u') and construct the associated curve 

v 1 = 4(A 2 - 3. 92 )(A + 3ei)(A + 3e 2 )(A + 3e 3 ), (21) 

The Hermite polynomial 3^s(p(t), A), depending on the argument t + u/, [IE [57| associated to the Lame potential 
6p(t + a/), which is already normalized has the form 

J 3 (p(t + lo'),X) = A 2 - 3p(* + lo')X + 9p 2 (t + w') - ~g 2 . (22) 

Using explicit expression for Hermite polynomial (|22[) we obtain the following simple solutions for the system (|9|): 

<h=-6p(t + w'), g 2 --2 4 -3 3 -p(i + ^) 2 + 2 2 -3 2 ff2 , (23) 

where a,Q = 0, ai = 3 • 4 • 7 • g 2 , 04 = — 4 4 • 3 4 • g 2 #3. More general solution with ^ can be written using the shift 
(HU in (HI]). Then we have 

gi = -6p(i + w') - -a , (24) 
q 2 2 = -4 2 • 3 3 p(t + u' f - Hp( t + w ') ao _ ± a 2 + 36ff2; (25) 
where ai = 3 • 4 • 7 • <?2 — 2ag/5 and 

1 5 84 3 432 2 !728 2 nn „ „ 
° 4 = 3125°° ~ 125 a ° 52 ~ ~25 a ° 93 + ~5~ a ° 92 + 20736 S2ff3- 

After changing the variables 

a = 2 2 • 5 • A, v = 2 s • a 4 (26) 

we obtain the curve of genus 2 

v 2 = 4(A 2 - 3 ff2 )(A 3 - 9A 52 + 27 ff3 ), (27) 
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where 32,33 are elliptic invariants (see for example 45]). The last expression is the another form of (f2Tj) . were 
Bj,j — 1, 2, 3, e 3 < e 2 < e 1 are real roots of equation 4A 3 — 32A — 33. The Weierstrass function p = p(t + u') shifted 
by half period ui 1 is related to sn Jacobian elliptic function with modulus k by [58[ : 

p(t + uj'; 32, 33) = a 2 k 2 sn 2 (at, k) + e 3 , (28) 



where a = \fe\ — e%. Using wave height a and modulus k — \/(e2 — e^)/(ei — e^) we have the following relations (see 
for example [45|): 

e 1 = 2-k 2 , e 2 = 2k 2 -l, e 3 = -(l + fc 2 ), 

32 = 2(e 2 + e 2 + e 2 ) = 12(1 - k 2 + fc 4 ), 

33 = 4 ei e 2 e3 = 4(fc 2 + 1)(2 - fc 2 )(l - 2fc 2 ). (29) 

From Lax representation for polynomial F we have the following nonlinear differential equation with spectral parameter 
A 

\FF - \P 2 (u(t) + X)F 2 + \ V 2 {\) = 0, (30) 

with eigenvalue equations 

v 2 (A) = 4A 5 - 21 A 3 9 2 + 27A 3 2 + 27A 2 53 - 8I3233 = 0. (31) 

VI. (n + l)-DEGREES-OF-FREEDOM GENERALIZED HENON-HEILES- TYPE SYSTEM 

We consider a generalized Henon-Heiles-type system with n + 1 degrees of freedom (2(| [H, [H, [33[ : 

1 - 

90 + 3g 2 + 2 E $ + = °> ( 32 ) 
i=i 

C 2 

q) + qaqj - -4 - a^j = 0. (33) 
"3 

where j = 1, . . . , n. Its Hamiltonian is 



1 " 1 ™ 1 1 ™ / 

# = 2 (p o + + «o + E 9? + 2 ao<Z ° 2 ~ 2 



2 



- ^ ) , (34) 



M 3^3 „2 I ' 

i=l 3 ~ 3=1 \ V I 

where Qo, 9j,Po,Pj, j — 1, ■ ■ - ,n, are the canonical coordinates and momenta, respectively, and ao, C 2 , a^, j = 1, . . . , n, 
are free constant parameters. The function H for n = 1 is the Hamiltonian of a classical integrable Henon-Heiles 
system with the additional term C 2 /g 2 . 

Next we will present (2 x 2) matrix Lax representation for the generalized Henon-Heiles system (|3~4"]) . The Lax 
representation has the form 

L=[M(X),L(X)}, L=^ U v y m =(qI) (35) 



where U,W,Q are 32, 33, 35]: 



(1 1 1 n 2 \ 

A+ 2 9o + 4 ao_ 16^1^) 

W = ~F + QF = a(A) U 2 - ^q X + Ja A 



-a(A) 
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j=i 13 

Q(t,X) =X-qo- 

The corresponding algebraic curve det(L(t, A) — |7) = of genus n + 1 is 

^ 2 = 4a(A) 2 ( A 3 + -a \ 2 + -a 2 A + -H + — V + — V , 

V M 2 16 8 16 ^A-ai 256 ^ A - a t 2 

\ 2=1 Z— 1 V 



(36) 



where 



,2 



1 / 2 1 1 1 2 2 \ 2 1 / 9 e 

Hi=- -PoQiPi - I a % - -a t q + -a a 4 + -q + ^ 1i I ft + ~qo \p\ 

e 2 9fe e 2 ft 



d2 



-a p, + aiPi + ya t + -a j 

7^ — - — ( (qiPk - QkPi) 2 - 



n 

a(X) =]J{X- a,), n>2. 

i=l 

Further we will collect the well known facts from finite zone inverse scattering transform method [3, ElJ 5H, H3, 
EH; H2> [fill useful for our construction of exact solutions of generalized Henon-Heiles-type system. Here we follow the 
Krichever construction [62|. The Baker- Akhiezer (BA)-function ip(t, A) of a nonspecial divisor D of degree n + 1 is 
given explicitly by 

A) = C(P) exp(iO lt ) Z) , (37) 

where #(u|.B) is a Riemann theta function, B = (B)ij = J b ujj is a Riemann matrix, u>x,u>2, ■ ■ ■ are normalized 

differentials of first kind, A(P) is Abel map Ak(P) — J P Wfc, <Zo is arbitrary point on real hyperelliptic Riemann 
surface Q36p denoted by K, Qi is the normalized differential of second kind with main parts at 00, dk, Z — tu + Z , 
Zo — —A(D) — Kq, Kq is the so called Riemann constant vector, the factor C(P) gives the normalization in (|3"T|) . 
The BA-function of nonspecial divisor D is the solution of 

j^-u(i))v&A)=A^(t,A), (38) 

where 

d 2 

u(t) = 2— logfl(Z) + const. (39) 
ere 2 

By the Riemann-Roch theorem there exists an nonspecial divisor D T and unique abelian differential f2 such that 
(D) = D + D T — 2oo and fi = (1 + 0(k~ 2 ) dk at 00. There exists a unique function ■0 T (£, A) called dual BA-fmction 
of D T . The BA-function ip T (t, A) of a nonspecial divisor D T of degree n + 1 is given explicitly by 

^(i, A) = C T {P) exp(-in lt ) 9{A f { ] z ~ Z) , (40) 

The factor C T (P) is fixed by the normalization. 

Let u(t) be a real finite-gap potential such that the Hill's operator has only finite number n eigenfunctions defined 

by 

tpi = atiip(t,pi), ipi = /?. i ^ r (t,p i )> 

atpi = Res Q.E, i = l,...,n (41) 

p=Pi 
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then the following expansion of u(t) in terms of squared eigenfunctions have the form 

n 

u(t) = - ai (3i tp (t, Pl )iP T (t, Pi ) + const. (42) 
i=i 

Let us construct the meromorphic differential 

f2 = Eip{t,P)ij) r {t,P). (43) 

A straightforward calculations give us 

Res fi = -u(t), 

P=Po 

Res n = -ip{t, P i)ip T {t,pi)a,l3 u (44) 

P—Pi 

where OLifii = Res ClE and Res Cl = — Res Cl. 

P—Pi P—PO P—Pi 

In particular case when K is hyperelliptic Riemann surface (j3"B"|) or in another form v 2 = 4]X™^~ (A — Aj) = P(A). 
The points of K are pairs P = (A, R) and A(P) is the value of the natural projection P — ► A(P) of if to the complex 
projective line CP 1 . 

For given nonspecial divisor D, there is an unique Baker- Akhiezer (BA) function ^(t, A), such that 

(i) the divisor of the poles of 'J is D, 

(ii) '5 is meromorphic on K\oo 
(hi) when P — ■+ oo 

oo 

*(t,P)exp(-fct) = 1 + Y J m s(t)k~ s , (45) 



is holomorphic and k = y/ A(P) is a local parameter near P = oo. 
There is a unique function u(x) such that 

* - u(t)# = A(P)f, (46) 
where \Jj is a BA function. Inserting expansion (|45[) into (|46[) . we obtain 

* - 2m(*)tf - A(P)# = exp(fet)0(fc -1 ), (47) 

and due to the uniqueness of 'J, we prove (|46[) . with m(x) = 2ro(i). 

By the Riemann-Roch theorem, there exists a unique differential Q and a nonspecial divisor D T of degree ri such 
that the zeros of f2 are D + D T and the expansion at P = oo, tl(P) = (1 + 0(k~ 2 ))dk. 

For given nonspecial divisor D T , there exists a unique dual Baker- Akhiezer (BA) function such that 

(i) the divisor of the poles of W is D T , 

(ii) '5 is meromorphic on fsT\oo 
(hi) when P — * oo 

oo 

* T (t, P) exp(-fci) = 1 + 51 ( 48 ) 

s=l 

Fix r to be the hyperelliptic involution P = (A, P) — * (A, -P), then we have D T = tD, ^ r (x, P) = *(x,rP). Let 
2<Li A*i(0) t> e the A-projection of D, and X^ILi MiW be the A-projection of the zero divisor of ^(x, P). The function 
^(i, P)^ 1 "^, P) is meromorphic on CP 1 and the following identity takes place 

$(l,Wt,P) = ?P (49) 

where P(i, A) = 117=1 — an d Mi W satisfies the following system of differential equations (Kovalevski-Dubrovin 

equations [6l|) 



4^W = 2 /^ -^y (50) 
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with initial conditions 

Mj(0) e [A 2i -i,A 2i ]. (51) 

Equations (|5D|) are first written by Sonya Kovalevski for genus n = 2 in relation to the integrable case of the Kovalevski 
top and for Korteweg-de Vries hierarchy of equations by Dubrovin for general n ( (6"H and references therein). These 
equations are useful for numerical calculation of the polynomial F(t, A). 
Introduce the Wronskian 

P), * T (t, P)} = P)t> T {t, P) - (t, P)* T (i, P) 



2 ^ (52) 



nr=i(A-^(o))' 

and the differential f2 is given explicitly by 

We assume that E(P) is a meromorphic function on i(T with n + 1 simple poles oo,pi, . . . ,p n and at P — » oo, E{P) = 
k + . . ., and P(P) is meromorphic function with n + 1 simple poles qo, <Zi, • ■ ■ , In and at P — > oo, P(P) = fc _1 + . . .. 
We also suppose that the divisors of poles of E(P) and E(P) are different from D, D T . 

^(t, A) = VfP exp jf ^^y^') • (54) 
A brief computation reveals that the BA function solves Hill's equation 

r 

and for dual BA function i\) T we have 



iH , u(t) ) ij = AV>, (55) 



r(t,\) = V / P(M) exp jf* ^r^y^' 



The solutions of the system with Hamiltonian in terms of Novikov polynomials F(t, A) of degree n + 1 in spectral 
parameter A 60(, for special points di i = 1, . . . ,n in closed intervals [A2«— ij A21], i = 1, . . . ,n the functions are given 
by miElll 

q = -u(t), g? = 16 / tlfl| — r, i = l,...,n. (56) 

u(t) is the famous Its-Matveev formulae ([59")) , [0], [63| and the points ttj lie in the lacunae [A2i-i, \2i\1 i — 1, ...ti, 
for generalized multidimensional Henon-Heiles system and are branch points in the case of the multidimensional 
Henon-Heiles system (p4|) with & = 0,j = 1, . . . ,n. The solution is real under the choice of the arbitrary constants 
di,i = 1, . . . , n in such a way, that the constants aj, i = 1, . . . , n lie in different lacunae of Riemann surface -fC . Then 
the constants Cj are given as 

e? = const 



2 ' 

a fc)! 



where i = 1, . . . , n and v is the coordinate of the curve (|36p .The constants are fixed by the initial conditions. 

Suppose that F(t, A) = ^(i, A), where 7 is Hermite polynomial associated with Lame potential u = (n + l)(n + 
2)p(i + lo'). Then the finite and real solution of the system ([52")l . ([55]) is given by (|5"1))) with the Hermite polynomial 
depending on the argument t + lo' (the shift in lo' provides the holomorphity of the solution). Consider the potential 
12p(t) and construct the associated curve [E3| 



v 2 



4A[|(A 2 ~6 ei A + 45e2-15 52 ), (57) 



10 



with branch points given by 

A = 0, Ai, 2 = 3 (l + k 2 ± 2^4 - 7fc 2 + 4fc 4 ) , 
A 3 , 4 = 3 f 1 - 2fc 2 ± 2^4~k 2 + k A \ , 



A 5 ,6 = 3 (fc 2 - 2 ± 2\/l ~ fc 2 +4A: 4 



(58) 



or in another form 



,2 



i/- , 7 63 - 297 , 4 4185 ,, 2 18225,, 
-j = x - y x 92 ~ T x 93 + ~le~ x 92 + ~ ~ X 9293 

/ 91125 2 3375 ,\ , 
+ {—9t ~ ^l) A- (59) 

The Hermite polynomial 3{p(t), A) associated with the Lame potential 12p(t) has the form 

2(p(t + u/), A) = A 3 - 6p(t + lu')X 2 - 3 • 5( 52 - 3p(t + u') 2 )\ 
3 2 • 5 2 

— (4p(iW) 3 -.g 2 p(t W)- 53 ). (60) 



Then the finite and real solution of the system (j32|) , (|33|) is given by (|56|) with the Hermite polynomial depending on 
the argument t + to'. The constants Ci, 62 are given by 

C 2 = const V{ai? 2 , C 2 = const.- (61) 
(ai - o 2 ) (a 2 - ai) 

where i = 1, . . . ,n and v is the coordinate of the curve (|5T|) . 

Next we list the periodic solutions of the system (f3"2"|) , (f3"3")l for n = 2 and for genus 3, where Ci = and 62 = 0. 
The (2n + 3) Lame polynomials of order n + 1 are solutions of 

^ + ((n + l)(n + 2)fc 2 dn 2 (ai) - A*) ^ = 0. (62) 

For n = 2 we introduce the following eigenfunctions E^ 3 -*, i = 1, ... 7 and eigenvalues Ai given in table HI Ai < A2 < 
. . . < A7 and the results are collected in Table HI see for example [5!| . 



TABLE I. 



1 


E< 3) 


= sn(crf, fc) 


Jdn 2 (^,fc) + C< 3) ) 


A< 3) 


= 7- 


5fc 2 - 


2^4- 


7fc 2 + 4k 4 


2 


Ef 


= cn(ai, fc) 


' dn 2 (at, k) +Cf ] ^ 


A< 3) 


= 7- 


2fc 2 - 


2V4- 


k 2 + fc 4 


3 


B?> 


= dn(crf, fc) 


[dn 2 (at,k)+Cf ) \ 


A< 3 > 


= 5(2 


-k 2 ) 




- fc 2 + 4fc 4 


4 


Ef 


= sn(ai, fe)cn(at, fc)dn(crf, fc) 


Af 


= 4(2 


~k 2 ) 






5 


Ef 


= sn(crf, fc) 


[dn 2 (at,k) + ci 3) ) 




= 7- 


5k 2 + 


2^4- 


7fc 2 + 4fc 4 


6 


Ef 


= cn(ai, fc) 


(dn 2 (at,fc)+Cf ) ) 


X (3) 


= 7- 


2k 2 + 2^4 - 


k 2 + fc 4 


7 


E 7 3) 


= dn(crf, fc) 


(dn 2 (crf,fc) +C* 7 3) ) 


A 7 3) 


= 5(2 


-k 2 ) 


+ 2vT 


- fc 2 + 4fc 4 



For convenience we present solutions in the following form 

90 = 12Q,fc 2 dn 2 (ai) ; g« = C®E® , = c£>Ef, ^j = l,...7 (63) 

where the constants Co, , C 2 are fixed by initial conditions. The same procedure is possible for general n. 

The Hermite polynomial 3 r (p(t), A) associated with the Lame potentials can be written in a different form, useful 
for applications: 

n 

■J(p(;),A)=^A fc (A)p(i)"- fc . (64) 

k=0 
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TABLE II. 



1 


C\°' = i(2k A - 3 - V4 - 7k 2 + 4k 4 ) 


2 


Cf> = i (k 2 - 3 - V4 - k 2 + k 4 ) 


3 


Cf 5 = i(2fe 2 - 4 - VI - fc 2 + 4fc 4 ) 


4 




5 


Cf > = ±(2fe 2 - 3 + V4 - 7k 2 + 4k 4 ) 


6 


Cf> = |(fe 2 - 3 + %/4 - fc 2 + k 4 ) 


7 


Cf> = i(2fe 2 - 4 + VI - fc 2 + 4k 4 ) 



For example for the genus 4 Lame potential 20p(i) we have 



2 6615 



A, = 11025, Ai=-1575A, A 2 = 135A^ —g 2 

A 3 = -10 A 3 + ^A 32 - 2450 ff3 

- 113, 2 3969 2 1925, 
A 4 = \ — — \ 2 g 2 + —gl + —*93- 

or in explicit form we have the Hcrmitc polynomial 3{p{t + a/), A) associated to the Lame potential 20p(i + u>') can 
be written as 

J(p(t + uj'), A) = 11025p0 + u'f - 1575p(i + uj'f\ 
+ (135A 2 -^ 2 )p(tW) 2 



1 QC r 

+ (-10A 3 + — A ff2 - 2450. 93 )p(t + w') 

, 4 113, 2 3969 2 195 , 
+A 4 - — A 2 .g 2 + —gl + — Xg 3 . (65) 



For n = 3 and genus four, the Lame curve have the following form 

/ 3 \ 



v 1 



4 1[{X 2 + 10e,A - 35e 2 - 7<? 2 ) (A 3 - 52A<? 2 + 560<? 3 ). 



V=l 

or in another form convenient for practical use we have 



,., x6 63129 l5 2 518505 , 4 
■j = A 9 - — A 7 <? 2 + ^-5sA 6 + -^-A 5 5 | — <?2S 3 A 4 

563227 o 4549125 2 \ , 3 991515 2l2 
16 - 52 + 16 53 J + ^~^ X 

f ^^5 2 4 - 5 -^^929l) A - 972405 335 3 - 1500625 ff 3 . (66) 



VII. SUMMARY AND CONCLUSIONS 



We approached the most general class among the three classes of integrable Henon- Heiles-type systems identified 




We provided exact quasi-periodic solutions for this class of systems, expressing these orbits via Kleinian hyper- 
elliptic functions. 

Applying the spectral theory for the Schrodinger equation to elliptic potentials, we pointed out elliptic periodic 
solutions. 

To obtain these results, we resorted to various mathematical methods, characteristic to both physics and astronomy. 
We applied these methods and results to a generalized Henon-Heiles-type system with (n + 1) degrees of freedom, 
pointing out the exact solutions in this case. 
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By emphasizing the exact solutions for the integrable class of systems under consideration, our paper contributes 
to a better understanding of the generalized Henon- Heiles-type problem. 
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